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The green mineral dioptase CueSieOig-BHhO has been known since centuries and plays an im- 
portant role in esoteric doctrines. In particular, the green dioptase is supposed to grant the skill 
to speak with trees and to understand the language of birds. Armed with natural samples of diop- 
tase, we were able to unravel the magnetic nature of the mineral (presumably with hidden support 
from birds and trees) and show that strong quantum fluctuations can be realized in an essentially 
framework-type spin lattice of coupled chains, thus neither frustration nor low-dimensionality are 
prerequisites. We present a microscopic magnetic model for the green dioptase. Based on full- 
potential DFT calculations, we find two relevant couplings in this system: an antiferromagnetic 
coupling J c , forming spiral chains along the hexagonal c axis, and an inter-chain ferromagnetic 
coupling Jd within structural CU2O6 dimers. To refine the J c and Jd values and to confirm the 
proposed spin model, we perform quantum Monte-Carlo simulations for the dioptase spin lattice. 
The derived magnetic susceptibility, the magnetic ground state, and the sublattice magnetization 
are in remarkably good agreement with the experimental data. The refined model parameters are 
J c =78 K and Jd=— 37 K with Jd/ Jc — —0.5. Despite the apparent three-dimensional features of the 
spin lattice and the lack of frustration, strong quantum fluctuations in the system are evidenced by a 
broad maximum in the magnetic susceptibility, a reduced value of the Neel temperature Tjv — 15 K 
■C Jc, and a low value of the sublattice magnetization m=0.55 \ib- All these features should be 
ascribed to the low coordination number of 3 that outbalances the three-dimensional nature of the 
spin lattice. 

PACS numbers: 71.20.Ps, 75.10.Pq, 75.30.Cr, 75.30.Et, 91.60.Pn 



I. INTRODUCTION 

Since ancient times, emerald has been one of the most 
rare and treasured gemstones because of its bright and 
brilliant green color. However, by far not all gemstones 
that were collected as emeralds or varieties of it were 
indeed emeralds — many of them later appeared to be 
specimens of green dioptase (see Fig. [J). Nevertheless, 
this confusion contributed considerably to the assignment 
of many mysterious powers to this mineral, such as pro- 
viding beauty, wealth and creativity. Especially, esoteric 
doctrines credit dioptase with the ability to grant the 
skill to speak with trees and understand the language 
of birds. Stimulated by such a rather complex realm of 
concealed powers, we attempt in the present study to 
unravel the also controversially debated magnetic prop- 
erties based on the paradigm of quantum mechanics and 
modern electronic structure theory. 

Dioptase is a copper silicate mineral forming remark- 
ably large shiny green rhombohedral crystals. Scientif- 
ically, it was first described and named by Rene- Just 
Haiiy in the famous "Traite de Mineralogie" in 1801. 1 
Vauquelin, as reported in Ref.Q], found that dioptase was 
a copper mineral containing silicate and - erroneously - 
carbonate anions. Only later, pure dioptase samples were 
analyzed and recognized as copper silicate with crystal 
water £ 

Structure determination showed that hydrous dioptase 
with the chemical composition Cu6Si60i8-6H 2 is a cy- 
closilicate with 6-membered single rings of silica SieOis 



(compare Fig.[TJ left) crystallizing in the space group i?3 
(S. G. 148)i2ii These rings are interconnected by Cu 2+ 
ions with a characteristic local environment of elongated 
octahedra formed by oxygen atoms. The Cu-0 bond 
length in this 4+2 arrangement is quite typical with 
about 1.96 A for the distorted equatorial plane of the oc- 
tahedra and a bit larger than usual (between 2.6 and 2.75 
A) for the apical oxygens belonging to the crystal water. 
Due to the sharing of the octahedral O-OH2 edges, the 
magnetic Cu 2+ ions form helical chains around the 3-fold 
axis along c (see Fig. [T|). Thus, each Cu atom has two 
Cu neighbors along these chains and another Cu neigh- 
bor with which it forms an edge-shared CU2O6 dimer that 
connects two adjacent spiral chains (see Fig. [T]). 



The magnetic properties of green dioptase have been 
investigated in some experimental studies^ - — Although 
these studies yield quantitatively slightly varying results, 
likely also related to dioptase samples originating from 
different locations, they essentially converge in the de- 
scription of dioptase as an antiferromagnet with a rather 
low Neel temperature (Tjv ~ 15 K) compared to the an- 
tiferromagnetic Curie- Weiss temperature of about 70 K. 
The ordered magnetic moment (m = 0.55/x.b) is dras- 
tically reduced with respect to the saturation moment 
of 1 [Ib for Cu 2+ . Together with the broad maximum 
in the measured magnetic susceptibility, 9 this puts the 
compound in the family of spin 1/2 quantum magnets 
that can be described successfully in many cases by the 
isotropic Heisenberg model 



H = ^2 JijSiSj, (1) 

at least for the low lying spin excitations. Here, Jij 
represents the exchange interaction between spins located 
at the lattice sites i and j. 

Although this model looks deceivingly simple at the 
first glance, neither its solution for a seemingly ordinary 
situation nor the assignment of appropriate exchange in- 
tegrals Jij for a specific material are trivial in any way. 
It is obvious that the crystal structure of a compound 
is the key to understand its magnetic properties. On 
the other hand, an assignment of interaction parameters 
solely based on structural considerations can be com- 
pletely misleading like in the case of (VO)2P2O7' 10 ' n In 
recent years, even careful investigations based on accu- 
rate experimental data, but within a limited spectrum of 
methods, have suggested controversially discussed mag- 
netic models for several compounds. A prominent exam- 
ple for this problem are the two closely related spin-1/2 
J\- J% chain compounds Li2Cu02 and LiCri2 02, for which 
consensus about their location in the magnetic phase di- 
agram was established only recently^ - — 

Thus, to establish the appropriate magnetic model for 
a new, complex material, the application of independent 
methods seems of crucial importance. In particular, the 
search for the relevant sector in the phase diagram can 
largely benefit from a detailed microscopic analysis based 
on modern band structure theor y 17 i 19 ~ — in combination 
with numerical methods to solve subsequently the corre- 
sponding Heisenberg Hamiltonian, at least in an approx- 
imate way 22 ' 23 

In particular, for green dioptase a magnetic model 
with antiferromagnetic (AFM) nearest-neighbor (NN) 
coupling J c along the spiral chains (see Fig. [TJ mid- 
dle) and AFM coupling Jd within the structural CU2O6 
dimers was suggested 2 ^ on empirical grounds and evalu- 
ated using quantum Monte-Carlo simulations (QMC) to 
fit the experimental magnetic susceptibility. This study 
places the compound in proximity to a quantum criti- 
cal point due to a competition between chain-like order- 
ing along c and magnetic dimer formation caused by the 
AFM Jd- In contrast to Ref. [24J, the results of our micro- 
scopic study place the compound in a different region of 
the phase diagram and assign the strong quantum fluc- 
tuations and the related magnetic properties to the small 
effective coordination number of the magnetic Cu 2+ sites. 



II. METHODS 

Electronic structure calculations were performed us- 
ing the full potential non-orthogonal local-orbital min- 
imum basis scheme f plo9 . 00-33^ For the scalar rel- 
ativistic calculations within the local density approxi- 
mation (LDA) the exchange and correlation potential of 



Perdew and Wang was chosen. A well-converged fc-mesh 
of 8x8x8 points was used for LDA calculations. 25 Wan- 
nier functions (WF) were calculated for the antibond- 
ing Cu 3d x 2_ y 2 states. Strong correlations are treated 
in a mean-field way within the LSDA+C/ approach. 27 
For the double counting correction (DCC) we applied 
the two limiting cases: the around-mean-field (AMF) 
approach and the fully- localized limit (FLL)j2£ The on- 
site Coulomb repulsion U^d was varied within the phys- 
ically reasonable ranges: [/3d=5. 5-7.5 eV for AMF and 
C^3d=6.5-9.5 eV for FLL. The intra-atomic Hund's cou- 
pling J%d was fixed to 1 eV. To allow for various spin 
ordering arrangements, the original hexagonal symmetry 
was reduced to the spacegroup PI. For the LSDA+J7 
calculations, we used k- meshes of 4x4x4 points. The 
calculations were carefully checked for convergence. For 
the structural input, we used the crystal structure from 
Ref. i. 

Quantum Monte-Carlo (QMC) simulations were per- 
formed using the programs looper and dirloop_sse of 
the software package ALPS.— The magnetic susceptibil- 
ity was simulated for ^=10752 sites clusters, containing 
256 coupled chains of 42 sites each. In the temperature 
range T/J c — 0.15-4.50, we used 25 000 sweeps for ther- 
malization and 300 000 sweeps after thermalization. The 
resulting statistical errors (<0.1%) are far below the ex- 
perimental inaccuracy. To evaluate the dependence of the 
static structure factor on the cluster size, we performed 
a series of simulations starting with a A=24 sites cluster 
and consequently increasing it up to 7V=8232 sites. Mag- 
netization curves were simulated on A=1536 sites clus- 
ters at T=0.025 J c using 50 000 sweeps for thermalization 
and 500 000 sweeps after thermalization. Statistical er- 
rors did not exceed 0.5%. 

The experimental data were collected on a natural 
sample of green dioptase. A green transparent crystal 
was mechanically detached from the calcite matrix and 
used for magnetic measurements without alignment in 
the magnetic field. The magnetic susceptibility was mea- 
sured with a Quantum Design MPMS SQUID in the tem- 
perature range 2-380 K in applied fields up to 5 T. 



III. RESULTS 

A. Electronic structure and magnetic model 

The electronic structure of the green dioptase was cal- 
culated within the LDA. The atom-resolved density of 
states (DOS) is depicted in Fig. H The width of the 
valence band is about 10 eV, similar to other cuprates. 
States at the Fermi level evidence a metallic solution in 
contrast to the green transparent crystals indicating an 
insulating behavior. This well known shortcoming of the 
LDA approach originates from the underestimation of the 
strong Coulomb repulsion in the Cu 3d shell. The insu- 
lating ground state can be restored by adding the missing 
part of correlation (i) via mapping onto a Hubbard model 
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FIG. 1. (Color online) Left: Crystal structure of the green dioptase. The CU2O6 dimers are shown in yellow and form a 3D 
network. SiC>4 tetrahedra are shown blue. The crystal water is depicted by O (blue) and H (gray) atoms with O-H bonds. 
Middle: the magnetic model of the green dioptase. Cu atoms are depicted as orange spheres, other atoms are not shown. The 
leading antiferromagnetic coupling J c (red) forms spiral chains running along c perpendicular to the projection plane. The 
ferromagnetic coupling Jd (blue) within the structural CU2O6 dimers couples the chains into a three-dimensional framework. 
Right: section of the spiral chain along c (bottom) and a natural sample of green dioptase grown on calcite (top). 



or (ii) in a mean-field way by LSDA+C/ calculations. In 
this study, both approaches are used. 

Despite the incorrect description of the ground state, 
LDA is known as a reliable tool for the evaluation of 
relevant orbitals and couplings. We start the analysis 
from the highest-lying states of the valence band. The 
well-separated band complex at the Fermi level is half 
filled and formed by antibonding Cu-0 dpa states. The 
energy range between —0.5 and —2 eV is dominated 
by non-bonding O and Cu states. At lower energies, 
around — 2eV, states of the SiC>4 tetrahedra and H2O 
appear. The sizable H2O contribution originates from 
the apical position of the crystal water in the distorted 
Cu04(H20)2 octahedra, and thus rather short Cu-Oh 2 o 
distances of 2.51 and 2.66 A. 

Typical for cuprates, the magnetic properties of the 
green dioptase are ruled by the half-filled antibond- 
ing Cu-0 dpa band complex at the Fermi level. The 
width of this complex W can be used as a rough esti- 
mate for the leading couplings. Thus, W=0.8eV for the 
green dioptase is comparable to related systems such as 
Li 2 ZrCu04 (buckled edge-shared chains, 1^=0.5 eV)^ 
Cu 2 (P0 3 )2Crl2 (distorted dimers, W=leV)^- or kapel- 
lasite (kagome lattice of corner-shared plaquettes with 
a Cu-O-Cu bond angle of about 107°, PF=0.9eV) 22 , 
but it is strongly reduced compared to S^CuOa (chains 
of corner-shared plaquettes, W—2.5 eV)^ 2 - or SrCu02 
(zigzag chains of edge-shared plaquettes, W=2 eV) 33 . 
Based on such simplified comparative analysis, we can 
conjecture the leading couplings in dioptase to be of the 
order of 100 K. 

The orbitals which are relevant for the magnetism can 
be evaluated by a projection onto a set of local atomic or- 
bitals. For each plaquette, one of the Cu-0 bonds and a 
direction perpendicular to the plaquette are considered as 
local x and z coordinate axes, respectively. This way, we 
find that Cu states in the dpa band complex have nearly 
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FIG. 2. (Color online) Total and atom-resolved LDA density 
of states for the green dioptase. The antibonding Cu-O dpa 
states form a well separated band complex at the Fermi level. 



pure Cu 3d x 2_ y 2 character (see Fig. [3J right), although 
the plaquettes forming the structural CU2O6 dimers are 
considerably distorted. On contrary, the orbital-resolved 
DOS for the O states shows a mixture of 2p x ,y and 2p z 
states. This mixing is caused by a non-coplanar arrange- 
ment of the neighboring dimers, sharing a common O 
atomi 3 ^ Thus, although the O 2p z contributions are un- 
usually high and seemingly hint at sizable O 2p v contri- 
butions, the states around the Fermi energy are clearly 
dominated by Cu-0 dpa states. Since the number of 
bands forming the band complex coincides with the num- 
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FIG. 3. (Color online) Top: Comparison of the antibonding 
dpa bands from the LDA calculation and the tight-binding 
model (left) together with the orbital-resolved density of 
states (right). Bottom: Cu-centered Wannier functions su- 
perimposed upon the leading superexchange pathways. 



ber of plaquettes in the unit cell, magnetic properties of 
the compound can be described by an effective one-band 
TB model. 

In case of the green dioptase, the evaluation of its mag- 
netic model from simple geometric considerations based 
on the crystal structure only is difficult due to a complex 
3D coupling of the structural dimers. To develop a mag- 
netic model of the compound from microscopic grounds, 
the six bands forming the complex at the Fermi level were 
mapped onto an effective one-band TB model parame- 
terized with transfer integrals tij . The Wannier function 
(WF) technique yields an unambiguous solution of this 
six-band problem. The resulting fit shows a nearly per- 
fect description of the LDA band structure (Fig. [3]). 

Only two of the resulting transfer integrals are rele- 
vant: t c = 126 meV (the subscript c stands for "chain"), 
running along the spirals of dimers (in the c direction), 
and td = 104 meV, the intra-dimer coupling (compare 
Fig. [1] middle). Other hoppings (except for the hop- 
ping ti C =24 meV between the spirals) are smaller than 
20 meV. To restore the insulating ground state, we map 
the TB model onto a Hubbard model considering an ef- 
fective (one band) Coulomb repulsion U c g — 4eV. For 
the strongly correlated limit at half-filling, both well jus- 
tified for the green dioptase, the lowest lying (magnetic) 
excitations can be efficiently described by a Heisenberg 
model. This way, the resulting magnetic exchange can 



be derived using the second order perturbation theory 
expression J AFM =4i 2 J /[/ c g . Since the original TB-model 
is a one-band model, only the antiferromagnetic contri- 
bution to the total magnetic exchange is accounted for 
in this approach. Thus, the resulting antiferromagnetic 
contributions for the leading couplings are J AFM =184 K 
and J^ FM =125 K. Since exchange integrals Jfj FM are 
proportional to tfj , all further exchanges are smaller than 
7 K (less than 4% of the leading exchange) and can be 
neglected in first place. 

Due to their close vicinity to 90° the intra-dimer Cu- 
O-Cu bond angle of 97.4° and inter-dimer angle of 107.6° 
require a careful estimation of the ferromagnetic con- 
tributions to the total exchange integrals, neglected in 
the effective one-band TB model approach presented 
above. Thus, we performed LSDA+C/ calculations of 
magnetic supercells with various collinear spin arrange- 
ments. The mapping of total energy differences onto a 
classical Heisenberg model results in an AFM exchange 
along the spirals chains J c =110 K and a FM intra-dimer 
exchange J d ——66 K for a typical value of f/ d =6.5 eV 
within the AMF DCC scheme. 

In agreement with the expectations according to 
the Goodenough-Kanamori- Anderson rules^~— the 
LSDA+J7 calculations evidence considerable FM con- 
tributions to both relevant exchange integrals J c and 
Jd- In particular, the total value of J c is strongly re- 
duced compared to the estimate from the one-band model 

(jAFM =lg4 R ^ yieldmg JFM = _ 74 K _ For j^^ the doser 

proximity of the Cu-O-Cu angle to 90° leads to an even 
larger FM contribution Jj M —~ 191 K which exceeds the 
AFM part J AFM , resulting in a significantly FM total 
coupling Jd within the structural dimers. 

Since the choice of DCC is non-trivial and can have a 
large impact on the resulting exchange parameter :] 23 ' 38 , 
we compare the AMF results to FLL calculations. For 
the green dioptase, we find that both DCC schemes yield 
similar results (Fig. 0]). The only apparent difference is 
related to the values of U^d- for FLL, about 2eV larger 
Uzd values are required in order to obtain the same ex- 
change integrals as AMF. The FM nature of Jd remains 
stable in the whole range of U^d and for the two DCC 
schemes. 

Although the qualitative microscopic model is well jus- 
tified by varying the U^d parameter in a rather wide 
range, the strong dependence of the resulting exchange 
integrals on U^d impedes an accurate estimation of the 
absolute size and the ratio of the two couplings. In the 
next section, we refine the values of the exchange inte- 
grals by alternative methods. 

To summarize the microscopic analysis, we obtain a 
model with two leading interactions: an AFM J c running 
along the spiral (in the c direction) and an FM Jd inside 
structural CU2O6 dimers. We should note that a related 
magnetic model was proposed in Ref. [24|. It is based on 
the same relevant exchange interactions, but implies an 
AFM intra-dimer coupling Jd in contrast to the FM na- 
ture of this coupling in our model. Remarkably, a model 



4 



140 
S 120 
100 



1 

- o 

- 

_ 


□ 

o 

! 


1 

□ 

o 

, 


1 - 

AMF" 
□ FLL ~ 

□ 

1 P 


1 


□ 


1 


1 1 1 






□ 




o 






□ 

□ _ 




o 


p 


i 


1 




1 


1 1 1 

i , ^< 



1 1 I I I I I I I I a 

5 6 7 8 9 10 

U d (eV) 



pare this model to the experimental results and refine 
the values of the leading exchange integrals. Since the 
microscopic model is well justified qualitatively (relevant 
couplings and their sign), its internal parameters J c and 
Jd can be refined by varying them in a reasonable range 
and subsequently simulating the thermodynamical be- 
havior for a given Jd/J c ratio. 

A method for simulations should be certainly consis- 
tent with the spin model. Since the two relevant cou- 
plings in dioptase form a non-frustrated, formally three- 
dimensional spin lattice, QMC simulations are natural 
and probably the only feasible choice. Therefore, we per- 
form QMC simulations for the relevant parameter range 
— }-<Jd/Jc< — 0.2 of the J c -Jd model. 

Simulations of the Heisenberg Hamiltonian yield a re- 
duced magnetic susceptibility x* which is related to the 
experimentally measured x by the expression: 



FIG. 4. (Color online) Results of total energy LSDA+t/ calcu- 
lations: the leading exchange integrals (J c and Jd) and their 
ratio as a function of the on-site Coulomb repulsion Uzd for 
the around-mean field (AMF) and the fully localized limit 
(FLL) double counting correction schemes. 

very similar to ours has been used to describe neutron 
scattering data for the dehydrated, black species of diop- 
tase CugSigOis. This issue will be discussed in Sec. IIVI 

B. Experimental results and model simulations 

To challenge our model with respect to the experimen- 
tal data, we measured the magnetic susceptibility of the 
green dioptase (Fig. [5]) and used neutron diffraction re- 
sults from previous studies^ Our susceptibility data are 
in good agreement with the earlier report o 9 ' 24 and show a 
broad maximum at 50 K along with the magnetic order- 
ing anomaly, evidenced by a kink at Tjv ~ 15 K. We do 
not observe any appreciable field dependence up to 5 T. 
Above 200 K, the susceptibility can be fitted to the Curie- 
Weiss law, modified with an additional temperature- 
independent contribution xo, responsible for core dia- 
magnetism and van Vleck paramagnetism. The Curie- 
Weiss fit leads to Xo=~ 7.2 • 10~ 5 emu/mol, the effective 
magnetic moment fj, e g=1.99 hb (5=2.30), and the Weiss 
temperature 0=43 K. Note that the susceptibility below 
Tjv depends on the crystal orientation in the magnetic 
field and can only be treated within an anisotropic model. 
In contrast, the magnetic behavior in the paramagnetic 
regime (above Tjv) is isotropic: the susceptibility sim- 
ply scales due to different g- values along different crystal 
directions ^ Since our microscopic study yields the ex- 
change integrals of the isotropic (Heisenberg) Hamilto- 
nian, we restrict ourselves to the data above TV. 

DFT calculations evidence two relevant interactions 
— the AFM intra-chain exchange J c and the FM intra- 
dimer (inter-chain) exchange Jd- Below, we will com- 



where Na, ks and \ib stand for the Avogadro constant, 
the Boltzmann constant and the Bohr magneton, respec- 
tively, g is the Lande factor, Ci mp is the Curie constant 
to account for possible impurity and defect contributions 
and xo is a temperature independent term, similar to the 
Curie- Weiss fit. 

Although a fit to x(T) is commonly regarded as a sensi- 
tive probe for internal parameters of a magnetic model^S 
we find that the ratio Jd/J c can be varied in a rather 
wide range (—0.8...— 0.4) yielding a very good fit to the 
experimental data for the paramagnetic phased. To im- 
prove our refinement of Jd/J c , we have to address the 
magnetic ordering temperature Tjv, which can be traced 
by a clear kink in the simulated curves. The reference 
to Tjv yields Jd/J c close to —0.5. The respective fit is 
shown in Fig. O (top). The resulting J c =78 K agrees 
well with the DFT estimates: 110 K for c7 3d =6.5 eV 
within the AMF scheme (Jd,/J c =— 0.6) and even better 
with 85 K yielded by U3d=8.5 eV for the FLL scheme 
(J d /J c =-0.55). Moreover, g = 2.26 and x =-6.9-10- 5 
emu/mol are consistent with estimates from the Curie- 
Weiss fit (2.30 and — 7.2-10 -5 emu/mol, respectively). 

For a further test of our model, we will address its 
ground state properties. First, the propagation vec- 
tor q of the AFM ordered GS coincides with the ex- 
perimentally observed q=(0,0, §7r)^ in the whole range 
— ]-<Jdl 'T c < — 0.2. In this GS, the neighboring spins 
along the spiral chains (J c ) align antiferromagnctically, 
while the ordering within the edge-sharing dimers (Jd) is 
FM. This justifies the validity of our microscopic model, 
but does not allow for a more accurate refinement of the 
proposed value for the Jd/ J c ratio. For a further compar- 
ison, we use the sublattice magnetization (m) that has 
been previously estimated in neutron diffraction experi- 
ments and amounts to 0.55 hb^ 

Unfortunately, the theoretical estimation of m is not 
straightforward for two reasons. First, the simulations 
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FIG. 5. (Color online) Top: The quantum Monte-Carlo 
(QMC) fit to the experimental magnetic susceptibility (per 
mol Cu). Bottom: QMC results for the ordered magnetic 
(sublattice) moment as a function of the Jd/Jc ratio. The ex- 
perimental value from neutron diffraction^ with an error bar 
is depicted by a striped bar. 
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FIG. 6. (Color online) Simulated high field magnetization 
curves for three different values of the Jd/Jc ratio. The mag- 
netic field is scaled by adopting the values of J c and g from 
the fits to x(T). 



to facilitate a further refinement of the model parameters 
due to practical resolution limits at high magnetic fields. 
In addition, the predicted value of the saturation field re- 
mains challenging for present-day experimental facilities. 



do not yield the magnetic moment in the ordered state 
directly. Instead, it can be estimated based on the static 
structure factor or on spin correlations. Second, 3D cou- 
pling is crucial for the magnetic ordering of the green 
dioptase, thus sizable finite size effects are expected even 
for rather large clusters. To account for these effects, we 
use the general procedure from Ref. and estimate the 
magnetic moment m based on a finite-size scaling of the 
static structure factor, taken for the propagation vector 
of the ordered structure. The results of simulations for 
various Jd/J c ratios are shown in Fig [5] (bottom). Re- 
markably, the theoretical m for Jd/ J c =— 0.5 is in good 
agreement with the experimental value. 

Finally, we can introduce a magnetic field term to 
our Hamiltonian and simulate the behaviour of magne- 
tization M as a function of the reduced magnetic field 
0<h*<5Jd- Such a simulation could be an additional 
test for our model, since high-field magnetization experi- 
ments were recently announced^ Therefore, we simulate 
M(h*) curves for J d /J c —~ 0.5 ±0.2 and scale them using 
the expression 

M(h) = M( f ^h*), (3) 
V 9VB J 

adopting the J c and g values from the fits to x(T). 
The resulting curves shown in Fig. [5] have similar shape 
and only slightly different values of the saturation field42 
Therefore, the experimental M (h) dependence is unlikely 



IV. DISCUSSION 

In our model, the spin lattice of the green dioptase 
comprises AFM couplings J c between the corner-sharing 
CUO4 plaquettes and FM couplings Jd between the edge- 
sharing plaquettes (Fig. [TJ. This situation is not sur- 
prising, because the corner-sharing connection normally 
leads to 180° superexchange, while the edge-sharing con- 
nection corresponds to the Cu-O-Cu angle close to 90°. 
However, the crystal structure of the green dioptase 
shows a tiny difference between the superexchange path- 
ways. The twisted configuration of the corner-sharing 
plaquettes leads to the Cu-O-Cu angle of 107.6° for J c 
that is substantially larger than 97.4° for Jd- The smaller- 
angle for Jd still fits to the general trend, predicted 
by Goodenough-Kanamori-Anderson rulesi2£r— On the 
other hand, the green dioptase is very close to the "criti- 
cal regime" of the Cu-O-Cu superexchange. Then, even 
a weak structural change could lead to a strong modifica- 
tion of the exchange couplings, making an empirical as- 
signment of the parameter region difficult. For example, 
the earlier theoretical analysis assumed AFM coupling 
for both J c and Jd^ Therefore, we used a quantitative 
microscopic approach and demonstrated that this em- 
pirical assumption is not consistent with the electronic 
structure of the compound. 

To examine whether small changes in the crystal struc- 
ture may lead to a modification of the microscopic 
model, the consideration of structurally closely related 
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compounds is a natural approach. The dehydration 
transforms the green dioptase into the black dioptase 
CugSigOi8 that essentially keeps the 3D framework- type 
crystal structure (Fig. [IJ but lacks water molecules. The 
Cu-O-Cu angles amount to 110.7° and 97.3° for J c and 
Jd, respectively— Thus, based on basic structural ar- 
guments, the signs of the two couplings should persist, 
while the absolute values are likely increased. This pre- 
diction is in line with the experimental data, indicating a 
large Weiss temperature 9 of 180 K 7 of the black dioptase 
compared to 0=43 K in the green dioptase (Sec. [Ill Bp . 
In addition, neutron diffraction studies evidence similar 
magnetic structures for the black and green dioptase. 7 
On the other hand, a very recent study of black dioptase 
based on extended Hiickel calculations and a ID fit to 
the magnetic susceptibility assigns the compound to the 
model family of AFM chain compounds with very weak 
interchain interactions. 45 

Further examples of the dioptase structure are 
given by the hydrated and anhydrous Cu germanates 
CugGegOig • iH^O with x — and 6. These compounds 
were previously considered as coupled frustrated spin 
chains, because a sizable next-nearest-neighbor coupling 
along the spiral chains was assumedi 46 ' 47 This assump- 
tion is rather empirical and mainly motivated by the 
chemical similarity to the well-known spin-Peierls com- 
pound CuGe03 with its frustrated spin chains of edge- 
sharing CuC>4 plaquettes4£ However, the pronounced dif- 
ference in the crystal structures strongly impedes a re- 
liable transfer of the well established magnetic model 
of the chain compound CuGe03 to the Ge-dioptasc 
CugGegOis; 46 Based on the results for the Si-dioptase, 
we would expect sizable AFM J c , while Jd is either FM 
or AFM. In the case that the FM and AFM contribu- 
tions to Jd are close to cancel each other, the inter-chain 
coupling is effectively switched off, and long-range cou- 
plings along the spiral chains could alter the physics. 
The above empirical analysis gives no clear evidence for 
significant long-range couplings, but further microscopic 
studies should challenge this conclusion. Thus, since the 
minor structural changes in the dioptase family might 
be crucial for changes in the leading magnetic couplings 
and the understanding of their magnetic properties, a 
detailed comparative study is underway- 
Taking the green dioptase as an example, we have de- 
rived the basic features of the dioptase spin lattice. This 
spin lattice is unfrustrated, hence we should preclude any 
references to the frustrated spin chain model, at least 
for the green dioptase CugSigOis • 6H2O. It is worth to 
mention that the dioptase structure does not give rise to 
the star lattice (decorated honeycomb lattice), as it may 
seem on the first glanced Such confusion could arise 
from a specific projection of the crystal structure, where 
the spiral chains look like flat frustrated triangles (com- 
pare to the middle panel of Fig. Q}. 

After shortly outlining what the dioptase spin lattice is 
not, it is more important to establish what it actually is: 
uniform AFM spin chains aligned along the c direction 



are arranged on the honeycomb lattice, i.e., each chain 
is coupled to three neighboring chains, and the system 
is geometrically 3D (Fig. [IJ. However, the total coordi- 
nation number is as low as three: each atom has two J c 
bonds and one Jd bond only. Thus, the couplings in the 
ab plane form a kind of a "sparse" honeycomb lattice. 
The reduction in the coordination number has strong ef- 
fect on the magnetic properties. 

Experimental data for the green dioptase evidence 
strong quantum fluctuations: the broad susceptibility 
maximum at T* ax / J c ~ 0.64, the low Neel tempera- 
ture {Tm/Jc — 0.2), and the reduced sublattice mag- 
netization (0.55 [Ib compared to 1 /xg for the classical 
spin-i systems). Strong quantum fluctuations are usu- 
ally observed in low-dimensional and/or frustrated spin 
systems. For example, the archetypal two-dimensional 
spin model of the square lattice reveals the susceptibility 
maximum at /J ~ 1.0 and a sublattice magneti- 
zation of 0.6 /ib-— To reduce the ordering temperature 
down to Tjv/J = 0.2, a very weak interlayer coupling 
J±/J ~ 10~ 4 is required^ Thus, the quantum fluctua- 
tions in the dioptase spin lattice are even stronger than 
in the square lattice, despite the 3D geometry. 

Quantum fluctuations in a 3D spin system can arise 
from the magnetic frustration (see Ref. 52 for an in- 
structive example). However, the dioptase spin lattice is 
neither low- dimensional, nor frustrated, hence its quan- 
tum behavior has a different origin. We suggest that the 
long-range magnetic ordering in dioptase is impeded by 
the low coordination number of the lattice, because the 
low number of bonds reduces the exchange energy that 
should stabilize the ordered ground state. The dioptase 
lattice can thus be compared to low-dimensional spin sys- 
tems with similar coordination numbers. For example, 
the honeycomb lattice having three bonds per site re- 
veals the low sublattice magnetization of 0.54 /j,b and 
T* ffl /J ~ 0.7 (compare to 0.6 /j,b and 1.0 for the square 
lattice with four bonds per site)^ The apparent simi- 
larity between the dioptase and the honeycomb lattice 
clearly shows that the coordination number is the actual 
criterion of the "low-dimensionality" , as long as the mag- 
nitude of quantum fluctuations (the tendency towards 
the quantum behavior) is considered. Although the con- 
clusion is a natural consequence of simple energy con- 
siderations, this is often overlooked. While neither the 
dioptase crystal structure, nor its spin model look low- 
dimensional, the essential physics is governed by strong 
quantum fluctuations, typical for low-dimensional mag- 
nets. The above considerations should stimulate further 
studies of dioptase-structure materials and the respective 
spin model. 



V. SUMMARY AND OUTLOOK 

Based on density functional calculations, quantum 
Monte-Carlo simulations and magnetic measurements we 
have derived a new magnetic model for the natural min- 
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eral green dioptase Cu6Si60i8-6H 2 on a microscopic 
basis. We have shown that green dioptase can be well 
described by a quantum spin 1/2 Heisenberg model with 
essentially two relevant interactions: an NN AFM intra- 
chain coupling J c ~ 78 K within the spiral chains run- 
ning along the crystallographic c direction, and a NN FM 
intra-dimer (inter-chain) coupling Jd ~ — 37 K within the 
structural CU2O6 dimers. The calculated temperature 
dependence of the magnetic susceptibility, the magnetic 
ground state, the ordering temperature and the sublat- 
tice magnetization for the suggested model parameters 
are in very good agreement with the experimental data. 
From our results we conclude that the dioptase spin lat- 
tice is neither low-dimensional nor frustrated, but ex- 
hibits large quantum fluctuations due to a small effec- 
tive coordination number of its magnetic sites despite 
the three-dimensional lattice geometry. 

Our approach demonstrates the great potential of the 
combination of modern band structure methods and nu- 
merical simulations with magnetic measurements for a 



reliable modeling of the magnetic properties for complex 
compounds. An empirically based assignment of inter- 
action parameters for structurally complex systems can 
be easily misleading and restrict studies to inappropri- 
ate regions of the magnetic phase diagram^ Since mi- 
nor structural changes may cause drastic changes in the 
leading magnetic couplings, especially for couplings via 
bonds close to 90° relevant in the dioptase family, a de- 
tailed comparative study for the hydrous and anhydrous 
Si- and Ge-dioptase compounds is in progress^ 
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